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We consider the case of a deformable material interface between two immiscible moving media, 
both of them being magnetiable. The time dependence of the metric at the interface introduces 
a non linear term, proportional to the mean curvature, in the surface dynamical equations of 
mass momentum and angular momentum. We take into account the effects of that term also 
, in the singular magnetic and electric fields inside the interface which lead to the existence 

of currents and charges densities through the interface, from the derivation of the Maxwell 
equations inside both bulks and the interface. Also, we give the expression for the entropy 
production and of the different thermo-dynamical fluxes. Our results enlarge previous results from 
other theories where the specific role of the time dependent surface metric was insufficiently stressed. 
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I. INTRODUCTION 



A lot of work has been already done on the study of the dynamics of a interface between two continuous 
immiscible media (gaseous/liquid or liquid/liquid). From a macroscopic point of view, the thickness of the 
interface is generally assumed to be negligible, with respect to its extent. A certain number of methods 
deal with this problem. The study of the interface between two media would be quite simple if one could 
consider the interface as a mathematical boundary problem. In that very first approach, the interface is not 
yet considered as a "real" physical phase by itself since it doesn't have any intrinsic quantity such as mass 
density, current density, or momentum density. If the electo-magnetic field is taken into account, one have to 
c/} • add the usual boundary conditions for the normal and tangential components at the interface between both 
bulks like for a ferrofluid in the analysis of the Rosensweig instability [17| . 

Even if very useful, this first approach, hides the true physical nature of the interface. Already, for pure 
fluids, the Laplace law expresses the jump of momentum balance along the interface to be equal to the surface 
7-H ■ tension times the mean curvature. Surface tension is the free energy of the surface: any real interface has a 
physical character leading to a much more complex problem. For example, how to take into account precisely 
a surface mass density, a surface flow ?And in the case of magnetizable media how to allow for the existence 
of a singular current between the twophases. To introduce explicitly a finite thickness is not a trivial matter. 
It leads to very tedious calculations [J H, d, 0, H- Indeed, to use this, for the electric fluid, see [20l |. or for the 
magnetic fluid without electrical effect, see — [j| becomes very difficult, due to the expansion in series of the 
dimensionless thickness, as for example for a conducting fluid with fully electro-magnetic effects. We will not 
follow this development any further. Thus, we need another approach where the interfacial physical quantities 
are introduced and are singular in the following way: they exist inside the interface and are meaningful as far 
• • . as the interface exists. They have no existence normally to the surface. Now the interface is a two-dimensional 
phase described in terms of intrinsic physical quantities for which dynamical balance equations have to be 
established, as is already the case for the adjacent bulk phases. This macroscopic as well as mesoscopic 
approaches permit still to ignore explicitly, but up to a certain point pj, the thickness of the interface [8(, 
even though it is different from zero, from a microscopic point of view. Indeed, we will suppose that any 
characteristic length is much larger than the characteristic width of the interface [7|, [9|, l20( 

A first entirely macroscopic approach was developed systematically in pjj [U, H3] , in terms of continuous fluid 
mechanics. The dynamical equations (mass conservation, energy conservation and momentum conservation) 
for the interface are obtained, using a 2 D physical model of the interface corresponding to the 3 D one de- 
scribing the bulk phase. The method consists in choosing a sample volume and in calculating for each physical 
entity its balance across the volumic sample, where the interface introduces a discontinuity whenever the 
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integration volume contains part of the interface. The integration introduces fluxes of any physical quantities 
which have to be exchanged between each bulk at their common interface. This procedure leads to two sets 
of balance equations: one for the bulk phases which is of course the usual one describing a single volumic 
phase and another one describing the interface itself. But that last takes into account the adjacent bulks since 
they are responsible for the fluxes crossing the surface. Those surface balance equations introduce a new term 
coming from the time evolution of the geometry of the deformable interface. This term appears already in 
the Laplace law, through the mean curvature. But in other surface balance equations, a term proportional to 
the mean curvature times the flux of the surface physical quantity along the surface appears. It reflects the 
local change in the metric of the surface itself and it is a highly non linear term which does not intervene if 
the interface is a pure discontinuity without physical characteristics. This fully "integral" approach has been 
used to establish the dynamics of the mass surfactant on a interface layer, see [24j and [26( and also for multi 
component system cf . [23j| . 

Another mathematical derivation has been proposed by D. Bedeaux, A.M. Albano and P. A. Wolff [E [Tol |. 
and was extended to polarizable cases by Albano and coworkers p], EE EI EE EH 

. Their aim was to build 

a suitable theoretical framework to apply non-equilibrium thermodynamics in fluid systems containing an 
interface and to relate it to statistical mechanics (TBI • In the mathematical basis has been done and they 
apply the theory to a two fluids system with interface where singular energy density is taken into account; the 
extension to the multi-component system was done by Kovac, [1 01 ] who also considered the case of fluid with 
internal angular momentum, (llj . However, this last author did not consider his model in conjunction with 
an imposed electromagnetic field. This author furthermore restricted his analysis to the flat undeformablc 
interface [lfj, since then the surface tension is independent from the choice of the dividing surface when one 
supposes local equilibrium [i"2| . 

The theoretical approach of Bedeaux and co. [E IE EE EH) EE EE EB| has been used for a system with 
electromagnetic effect: firstly without singular magnetic field in the interface, (llj and with singular densities 
and current in , EB]. But, each time, the balance of angular momentum has not been considered. The 

importance of the approach is how to treat the thickness: the interface position is considered as proportional to 
a delta function and the two bulk phases existence is taken into account formally by a Heaviside function, this 
formal decomposition has the advantage to suppress the problem related to the integration of the thickness. 
In one of their first papers, Albano [l4j discussed the very sensitive problem linked to the interface mean 
position Q and gave several reason to justify this method for a lot of systems. The change of shape of the 
surface is expressed by the mean curvature Ti. which appears explicitly only in the normal component of the 
jump of momentum at the interface. The change of metric expresses the local variation of distances along the 
interface: expansion or contraction of the surface. It has to be taken into account, as well as local stretching 



241. The non linear term related to the time dependent metric has been shown to exist in previous works 
19 , HE HE HE HE H3 • in fact, it appears explicitly in the equations, but it is not sufficiently stressed in the 
original approach of Bedeaux and coworkers, leading perhaps to some misunderstandings. In conclusion, all 
three approaches 0, EEEE HE HI HE HE Hi or Gogosov and al. [H, 0, S, i, H and the one derived by Bedeaux 



et al. [11|, |13|, [14|, Ha, [la] lead to the same surface balance equations 



This paper is organized as follows: in section two, we recall the essential mathematical basis of the theory, 
since much of the detailed derivations are to be found in Q . In section three, the dynamical equations for the 
linear momentum and for the electro-magnetic field, in the interface and for the bulks phases are established. 



Indeed, our main aim is to extent the basic study of Rosensweig 27J to the case where a free magnetisable 
material interface separates two immiscible fluids, neglecting relativistic terms [IE In section four we 
compute the governing equations for the total, internal and external angular momentum. In section five, we 
derive the constitutive equations for the magnetization and for the polarization. In section six, we establish 
the conservation equations inside the interface and the bulk for the kinetic energy, the rotational kinetic 
energy and for the electro magnetic energy. In section seven, the corresponding equations for the internal 
total energy are obtained. However, in each case, we sill stress explicitly the role of the time dependent metric 
since Bedeaux et al.[9] did not. In the next section, we discuss the entropy production and associated fluxes. 
We end up with the conclusion. 



II. MATHEMATICAL MODELING OF THE INTERFACE 



In this section, we present the mathematical basis to model the system, starting from the mathematical 
development explicitly done in the past[E IE EE E3|- We will review the general definitions of the different 
quantities to be used and then, we show how the time dependence of the curvature can be taken into account 
in this formalism, we will apply this new feature to the equation for the conservation of mass of the system. 



A. Mathematical Basis of the theory 



Let us recall the presentation made in [9J], also used by [10|. For simplicity, we consider two distinct and 
immiscible one-component fluids separated by a time dependent and spatially deformable interface, S(t), 
denned by the function f(x, t) = such that f(x, t) > corresponds to the fluid I and that f(x, t) < defines 
the fluid II. Each physical quantity, A is written as the sum of the two bulk phases quantities (A + for the 
phase I, A" for the phase II), and a singular quantity representing the interface itself, A s as: 



A(x, t) = A'®- + A + Q + + A S 5 S 



(1) 



Equation {T]) introduces two new functions which are the most important mathematical tools to be used. To 
begin, the Heaviside step function, defined as follows: 
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The Heaviside step function enjoys a very important property 



df 



±5(f) 



(2) 



where 5(f) is the Dirac function. Since exists a direct relation between the Heaviside function and the delta 
function, the "surface delta function", S 3 , which appears in fl} will be expressed in terms of 5(f ). To do that, 
Bedeaux et al. [9J] consider the total area of the interface: 



S(t)= / Vf(x,t) 5(f)dx 



The normal vector at the interface, n, pointing from bulk II into bulk I is defined as: 



V/ 



V/ 



(3) 



(4) 



From expression §S§ and ([¥]), Bedeaux et al. defined the function 5 s (f) as: 

6 s (x,t) = Vf(x,t) S(f(x,t)) 



(5) 



That function S B (x, t) has an enviable property: one can restrict any integration of the volume that contains 
the interface of discontinuity , to an integration over the discontinuity surface S(t). Using equations JHI31H]), 
we have: 



dQ ± (f(x,t)) 
dx 



±n6 s (x,t) 



By definition, the function / is equal to zero at the interface, so that we have: 

^ = -u s .Vf = -< |V/| , for / = 



(G) 



(7) 



where u s is the velocity of an object intrinsic to the surface and thus, it* = n.u s is the interface velocity, taken 
normally to the interface. It can differ from u il \f = Q.n which is the normal component of the ± bulk velocity u^ 1 
along the surface, like for example in non equilibrium crystal growth [32| . The difference [u±\ f{1 t) _ -n — u s n ].n 
defines thus a relative velocity 

a t tne interface and thus a flux towards or leaving it. 



The time derivative of the Heaviside function is: 

dG±(f(x,t)) 
dt 



Tu s J s (x, t) 



Like Bedeaux one defines the total time derivative at the interface f(x, t) — as follows: 

D s d 



Dt dt 



u s .\7 



(8) 



(9) 



From last equations Bedeaux et al. Q derive the important results: 

D s 9 ± (f) D S S S , , 

— ^ = — ^— = (10) 
Dt Dt y ' 

Another important property concerning the quantities on the interface is the following: 

n.VA s = (11) 

where A s is a physical quantity intrinsic to the interface. Equation (fTTjl means that the quantities of the 
type A a which are defined on the interface cannot change in the normal direction to that interface: they can 
change only along the interface. In practice, thus, whatever A s we consider, we have V„4 S as V S .4 S with the 
definition of the operator V s = [(I — nn).V]. 

The properties (|6|8|10|lip will be very useful in any future calculation. Indeed, the delta and Heaviside 
functions lead to a simple definition of the interface. To compute the equations for the bulks and the interface 
are simpler than in another approach for which the definition of the different phases in terms of boundaries 
in the integrals leads to very tedious calculations even for situations simpler than those that we are going to 
deal with [10, 0,1111. 



B. Influence of the time dependent metric: Example of the continuity equation for the density 

The influence of the time dependent metric is actually well known [Hj] — (23|. For example, in [24j and [2(|, a 
simple derivation of the dynamics of the mass surfactant at the interface under the influence of the motion 
of the interface is given. It reflects the time dependence of the compressibility of the interface. However, the 
works [H — take explicitly into account only the effects of the time depending mean curvature. The change 
of shape is considered by Bedeaux et al., but it lies hidden in the formalism. Indeed, as we stressed in the 
introduction, Bedeaux and coworkers 0, [13, [ill 0) EH take only into account the usual Laplace law, in the 
balance of momentum along the interface. In this section, we show how the time dependent surface metric 
can appear naturally, extending thus the mathematical formalism developed by Bedeaux and coworkers and 
finding back the results introduced by [ljj H3, [2l|, [H, [H, [24|, [26| , as part of our extension. 



As an illustration of this approach, we consider the total mass conservation law and derive its formulation for 
each bulk and inside the interface S. Let us consider the total density px, using the decomposition (fTJ), we 
have: 

p T (x,t) = p + (x,t)Q + + p-(x,t)e- + p s (x,t)S s (12) 

where p ± are the bulk densities of the ± bulk respectively extended up to the interface and p s might be the 
Gibbs superficial mass density, depending on the choice of the mean surface position. If we note the total 
volume of the system, Vr, we can rewrite it as Vr — S © V + © V~ . Supposing that the bulk phases are 
immiscible, we write the conservation of mass density in the integral form, in the absence of any chemical 
reaction: 

^ / p T dV T = (13) 



Dt _ 

If we put the decomposition (|12jl in (fT"3|) . we obtain 



^ J (p s 5°)dV T + jLJ (p±0±)dy T = (14) 

The first term in the LHS of (|14[) is proportional to 6 s . The next term is proportional to the Heaviside function, 
which is non zero only for the corresponding bulk phases, then we can rewrite (|14j) as 

^ J ( P s 8 s )ds + §-J ( P + e+)dv + + §- t J (p-e-)dv- = o (15) 

In the first term of (|15[) . the density p s is defined only on the interface, then we can write for this term 

Then the first term in eq. (fT5|) is 

^ J (p s 5°)dS = J %- t { P s 5 s )dS + J {p s 5 s )%- t dS (17) 



and now using the second relation in ea. (|10|) we have 

g J (p^)dS = J (^£)S s dS + J ( P s 5 s )^dS (18) 



In the RHS of (118|) . the last term is proportional to times 5 s . To evaluate it, we use the well known 



property [191, |2fl l2il, 122, |23| 

ddS 



dt 



J_dg 

2g dt 



dS so that 



6* A' 



ddS 
~~dt 



5 s A s — ^-dS 
2g dt 



(19) 



for any quantity ,4 s defined on the surface. By definition, g = det(gij), where gij is the metric on the interface 
and where [HM M, H, [H| : 



9 9 A -,B -fn. 



(20) 



In last expression ([20]) . we used the classical definition of the mean curvature H = -|V s .n (cf. [T||, Ql| [13, 

[11 11 11 [11 [II). 

Another way to obtain this important point rests on the demonstration given by Stone [24j ], who performed 
some mathematical transformations on the second term of the RHS of (fT5[) (see f [24l |25|) for more details): 



El 
Dt 



dS= [dSV}.u s - (Vu s ).dS 



where dS = ft.dS. Taking the scalar product of eq. (|2ip with n, one obtains 

D 



and because 



we finally obtain [251 ] : 



j^dS = dS(V.u s ) - (Vu s .n).ndS 



(Wu s .n).n = (nn.V). 



E^dS = dSV s .u s 
Dt 



Now, we can rewrite (pT8|) as follows: 

( P s S s )dS = 
we use the following result : 



Ei 

Dt 



dt 



u s .S/ s p s + p s V s .u s + p s V s {nn.u s 



dS 



p s V s (nn.u) = p s (V s .n){n.u) 



We finally obtain 



D 



^ / (S s p s )dS = 



= s 



+ u s .V s p s + p s V s .u s + p s (y s .n)(ri.u s ) 



dS 



and taking into account the above definition of the mean curvature Ti. , we have our main result: 

~dp s 



El 
Dt 



{S s p s )dS 



at 



u s .V s p s + p s \7 s .u s - 2TLp s {n.u s ) 



dS 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



The last term in the RHS of (|28[) introduces the mean curvature TL times the normal component of the surface 
velocity , enlarging the original demonstration given by Stone [13], since we are considering also the total 
mass flux coming or leaving the adjacent bulk phases. Furthermore, Stone did not give any clear definition of 
the velocity. This demonstration has been also used in the case of and Q. But, these authors considered 
that n.u s — it* the bulk velocity at the interface is the surface velocity. On the contrary, we will assume that 
the surface velocity at the interface is u B , but this does not mean that at the interface, the bulk velocity normal 



component u.ft calculated for f(x,t) = has the value ft.u s 



" drift" terms of each phase which are proportional to ft. [fp 1 1 



l/=0 



u s n of the surface velocity. In our approach the 
—u s \ are not neglected, as we will see just below. 



Now, we assume that each phase is incompressible. Thus, the total volume of that phase remains constant: 
0. The effects of the curvature are negligible inside the bulks and we have for the second term 



DV+ _ DV + 

n t Dt 

of (IM|) . using the Reynolds transport theorem (see [27|) and the relations ([6]) and 



D 
Dt 



( P ± e ± )dv T 



We can rewrite last equation (j2"§j) as: 

^ f { P ± e ± )dv T = fdv 



d P ± 
dt 



d P ± 



T ft.u s p ± 5 s + ^V. (iPp*) ± rt.u ± p ± S s 



+ 9 ± V.(u ± p ± ) 



J dS [Trt.u s p ± S s ± n.u ± p ± S s ] 



(29) 



(30) 



Thus summing (f2"5)) and ([3H|) . we obtain the explicit form of the integral equation (p~5j) . But, we know that 
the coefficients of 5 s and ± must vanish separately (cf.[^]), and thus we obtain the following results for the 
bulks phases: 



dp ± 
~~dt 



V.(5 ± p ± ) = 



and for the interface: 



DsP s 
Dt 



p s V s .u s - 2p s u s n H + p + ft.{u+ - u s ) 



.{u- -u s ) = 



(31) 



(32) 



This typifies thus the approach of Bedeaux and coworkers @, 0, [ljj, [HI, EH- Stone's physical explanation of 
the term proportional to — 2u^TL was thus extended to consider the influence of the drift terms, p + ft.(u + — 



u s ) — p~n.(u~ — u s ). In (|32|) . V s .u s — 2 u s n H expresses "a source like contribution resulting from local change 
in the area which are related to stretching and distortion" [24| . These two terms appeared already in (|25p since 
we are writing u s as a the sum of two contribution: one is along the surface and the other is normal to the 
surface. The first one shows up as V s .u s and expresses the compressibility of the surface along itself and the 
other —2 u s n T-L is due to the change of shape of the interface. The equation for the conservation of the mass 
inside the interface (|3"2"|) is exactly the expression derived in [24]],[26j or in the works by V.V. Gogosov and 
co-authors [lj, [2j, 13J, |4|, l5( if one supposes, as these last do, that the normal component of the bulk velocity at 
the interface is the velocity of the interface. 



III. DYNAMICAL EQUATIONS FOR THE INTERFACE 

In this section we consider the derivation of the dynamical equations of our system, using the method described 
above to show explicitly the effect of the curvature. We suppose that both bulk phases are pure liquids since 
the multicomponent bulks with chemical reactions has been the subject of extensive studies [ToL fl3l fl9L I20I l2ll. 
HH, [23j]. We will consider a system where magnetic and electric effects exist. Thus, we define the electric field, 
E, the magnetic field, H, the magnetic induction, B, with the relation: B = po(H + M), where M is the 
magnetization of the system and po the magnetic permeability; D is the displacement field which is related to 
the electric field through the relation D = EqE + P, where P is the electric polarization and Eq is the electric 
permitivity of the void. We define the charge density p e and the current density 7. However, for the time 
being, we will not consider terms linked to the semi relativistic approximation [29i . |30( and will neglect terms 
± k s 

proportional to — , — . The Maxwell equations are written in MKSA units, like for Rosensweig 1271 129| 
c c 

A. Momentum balance equations 

Consider now the momentum balance equations starting with the Navier Stokes equations, written in an 
integral form, see [U[H,[33|. The general expression of the momentum balance is: 



/ \ P u)dV T + J \p em )dV T = J[(V.f)+ P F ex ]dV T (33) 



where we have the total stress tensor T ( taking into account the pressure the viscous effects and an electro- 
magnetic part [13, Hil), the vector for electromagnetic impulsion p em = D x B and the external force of non 



electromagnetic origin (e.g.: gravity), F ex . Since we wish to extend the initial approach of Rosensweig [271 l29j| 
to the present case, we are using like him the Minkowsky expression for the electromagnetic momentum, while 
Albano and co [H, [lj, [l5|, [l6| , or Felderhof and Kroh [3(| took the Abraham form \ E x H in their derivation. 

Like we did for the density (fl2|) , we use the decomposition (fTJ) for the momentum, the tensor and the external 
force, replacing in (fl]) A{x, t) by pu, p e m, T, pF ex respectively so that (|33|) becomes: 



D 

151 



[p+u + e + + p~u 



■e-] dv T + ^- t J ( P s u s s s )ds + J [p+ n e + +p em e-] dv T + §- t J (p: m s s )ds 



(34) 



/V.(T 6-+T B+ + TS s )dV T + J(p-F- x O- + p+F+e+ + P s FiJ s )dV T 



1. Bulks and surface momentum 

To begin, let us calculate the three first terms appearing in the left hand side of last equation (j3"4f , that is we 
develop the time derivatives for the bulks and the surface momentum. We have 



D 

m 



[p + u+e+ + p-u-Q-} dV T 
dp + u + 



dp + u + Q + . I f t , t — ^ I . dp u 6 . _ _ _ , 

+ v.(p u u e> ) + hv.(^rre _ ) 



dt 



+ V.{p + u + u + ) 



Of 



e + dv + 



dp u 



Of 



dt 



e-dv 



+ j [p + u + n.(u + - u s ) - p u n.(u ~u s )]S s dS 

where we used the properties of the Heaviside and delta function (|8I9I10I11[) . 
The next term to be calculated at the interface is: 



D(p s u s 

dT 



dS+ / [p 8 u s 6 



d{p s u s ) 



Of 



+ V s .(p s u s u s ) ~ 2H(n.u s )p s u s 



DdS 
Dt 

6 s dS 



This result follows directly from the method defined in (|19[) . which was used to obtain 



(35) 



(36) 



2. Electromagnetic impulsion 
Let's consider now the term for the electromagnetic impulsion. Using (I19p . we obtain the following expression: 

^ / p^dVr + ^( p e s m S s dS = 



d(p e ±® M 
dt 



+ V.(u ± P e t® ± ) 



dV T 



d(S s p: m ) 



dt 



+ V.(u s p e s m S s ) 



dS 



2H(u s .n)p e s m 5 s dS 



(37) 



3. Stress tensors 

Now, we will calculate the RHS of eq. ([M|) Following |T]), we have to calculate the stress tensor: 



V.TdV T = 



v.(t e~) + v.(T e+) + v.(T(P 



dV T 



At this stage, it is useful to define a new stress tensor T 

T = f - Up e 



(38) 



(39) 



It is by definition, the total stress tensor minus the flow of the electromagnetic impulsion [27j , [29| . This enables 
us to obtain a much simpler expression for the momentum balance, since it eliminates term linked to V .up ern . 
Then, Eq.® gives : 



J (V.T + )0 + dV + + J (V.T )Q~dV~ + J V.T S +n.(T + -T ) 5 s dS + J(T.V6 s )dV T (40) 
Using {SJ and (fit))) one obtains the following relations for the electromagnetic impulsion term: 

^n^ = -P S e m {u s .V5 s ) and g" t = ~§f Q + + P™)*" (41) 

Inserting all the previous results ([35| . (|36| . (|37|) . (|40| . (|41|) in (|34|) . and taking into account that the term 
proportional in and S s must vanish separately, the momentum balance for each bulk is given by: 

^ + VV^) + °& = V.T* + P±Ft (42) 
while for the interface, its specific momentum balance reads 



(43) 



+ V.(pW) + ^ff - 2H(u s .n)(p s u s +jf em ) + (u s .n)(p- m - p+ m ) = 

n.{T + - f ~) + p~u~{u s - u s ).n - p+u + fi.{u + - u s ) + V .T'" + p s F° x 
We have also some terms proportional to V(S S , so that on the normal components: 

pl m u s .n + f S .n = (44) 



4- Explicit form of the stress tensor T 

The tensor T complete form, for each bulk has been introduced already by [13, [H, HO] : 

T = T + T em (45) 

Here, the tensor T is, as follows: 

f = T + 2£Z.(w-fi) (46) 

where £ is the vortex vorticity and Z is the alternating polyadic corresponding to l x l y l z eijk where 

(ijk = 123,231, or 312) 
(i = j,i = k, or j = k) 
1 (ij jfe = 132, 213, or 321) 

In Eq. PB|) we take into account of the effect of the vorticity (defined as uj — x u) as well as the internal 

angular momentum contribution Q. Those two last effects have not been considered in Q, or in (l3l. IT3. ll5L fl^| . 
Kovac, pH . considered the intrinsic angular momentum of individual particle rotation and thus spin, while 
ignoring possible electro-magnetic effects. In (l4l)l) . we will suppose that T has the following form describing a 
Newtonian fluid: 

T = -pi + 2r/D + A(V.u)I where I is the unit tensor and D = ^(Vu + (W) T ) (47) 




In the definition (|47j) . p is the hydrostatic pressure, 77 is the shear viscosity, A the volume viscosity. We assume 
that both bulk phases are incompressible, so that V.t? = 



We consider also the usual form of the Maxwell tensor [l 

1 



T 

1 em 2 



e E 2 + ti H 2 



l + DE + BH (48) 



As before, all of the physical quantities (velocity field, u, internal angular momentum,!},... ) are decomposed 
using the formula JT]). For the tensors, we take the following form: 

T em = T- m e- + T+ m e+ + T s em S s (49) 



with the definitions: 



rp± _ _1 

e m 2 



£ 0(S 2 )±+ M0 (ff 2 )± 



rpS _J_ 

em 2 



I + D ± E ± +B ± H ± 

(50) 

I + D S E S + B S H S 



We will define now a very simple form for the tensor T . Even if the adjacent bulks are incompressible, the 
material interface is compressible in the case where there is no drift. Thus, the description of the interface as 
a kind of 2D Newtonian liquid is somewhat more complicated than the one for a bulk fluid. Basing ourselves 
on previous results, [T3| the derivation leads to the following expression for the stress tensor at the interface: 



~8 



T = crl s + t s + 2f Z.(w" - ft s ) (51) 



with the intrinsic surface stress r s = 1rfT) s + A s (V s .u s )I s [I8l.ll9tl22l.l23l]. Let us note that the interface is a 
compressible 2D media if u s n — u.ri. Then we have 

V.f S = 2aHn + V s a + V s r s + V.(2£ s Z.(w s - &)) (52) 



The first term in the RHS of (j52|) is the Laplace force, proportional to the surface tension a times the mean 
curvature, Ji. This term is present in all previous works (cf. [I3|,[ll|). The second term introduces the variation 
of surface tension along the interface and expresses thus a Marangoni effect while leaving free its precise 
physica-chemical origin [l(| EH, E3] ■ We are not aware that the surface angular momentum (2£ s Z.(aJ s — f2 s )), 
introduced in the last term on the RHS of (|52|) . has ever been taken into account. Thus in a way, the total 
stress tensor has a very simple but heavy mathematical form 

B. Maxwell's equations in a moving medium 

We use the same kind of approach to derive now the Maxwell equations for the interface and for the bulks 
phases. Thus applying to the Maxwell equations the decomposition ([T]), we will obtain terms proportional to 
V(5 S , in addition to the ones proportional to <5 S and [1J] and [15fl. Like in previous derivations, this kind 
of terms will be related to the boundary conditions at the interface jsj . 

1. No explicit time derivative in the Maxwell equations 

For those two Maxwell equations, the derivation typified by {T} is identical to the ones given in [TJ] and [HI], 
since the time dependency of the surface metric (|19[) will then not intervene. We will show thus directly the 
results. 

A] Writing the equations in their local form the Gauss law: 

V.D = Pe (53) 

Using the decomposition (TT|) for the vector field D and the charge p e , and using the properties for the time 
and spatial derivatives of the 5 s and functions, we arrive to the following results: 

• terms proportional to ± : 

Eo (Vi ± ) + V,P ± =p± (54) 

• terms proportional to <5 S : 

e Q n.(E+ -E-)+ n.(P+ - P~) + V.(e E s + P s ) = p s e (55) 

• terms proportional to \7S S : 

n.(e a E s +P S ) = (56) 

B] The next equation is the zero divergence constraint on the magnetic field: 

V.B = (57) 
which gives us the following equations for the interface and the bulks: 



terms proportional to 0^ : 

V.B* = (58) 

terms proportional to 6 s : 

n.{B + ~ B-) + V.B s = (59) 

terms proportional to \7S S : 

n.B s = (60) 



2. The time derivative Maxwell equations 

Now, we will consider the time depending Maxwell equations, i.e. the Ampere law and the Faraday law. Since 
we assume that both bulks and the surface are moving at maybe different velocities, we consider that they 
are seen from a stationary fixed laboratory reference frame. We will calculate the local form of Maxwell's 
equation, starting from the total time derivative of the volume integrals or of the surface integrals and use 
thus in that process explicitly ([2"5"|) [HI, HI] . 



A] We begin with the Faraday law, whose integral form reads: 

D 



V x E.dS T = / B.dS T 

Dt 1 



(61) 



This expression is a relation between the curl of the electric field along a arbitrary surface St which encloses 
the volume sample ± and the interface between them, and the temporal derivative of the magnetic flux which 
goes through the surface St- As has become usual by now, we must take the total time derivative in the RHS 
of (|61|. since the surface S^deforms itself with time [H, [34]]. We use the same kind of formal decomposition 
for the dSr as we did for the dVr- 



dS T = dS + © dS~ © dS s 



(62) 



As we did already for the dV , we neglect the temporal variation of the dS surface elements, because both 
dS are chosen as arbitrary surface elements, inside the corresponding ± bulk. The situation is quite different 
for the interface dS s : for that last one, we must take into account its time variation. Starting from eq. (|61|) we 
have using the decomposition |T]) and the relation (|62| : 



V \ ( E (-) ' - £••<-> ' - E*,r ) .dS=-— 



Applying ©, © in this equation leads to 



B + Q + + B-Q- 



.dS- -jy I n-.r.dS" 



(63) 



J (V x E+)e + dS+ + J (V x E-)Q~dS- + J ((E~ - E+) x n + V x E s )S s dS s - J (E s x \75 s )dS 



D 
Dt 



B+0+ + B-Q- 



.dS- 



D S B S 
Dt 



S s dS s + {B s 5 



DdS s 
Dt 



(64) 



Now, let us look at the first term in the right side of equation . Following a procedure analog to the one 
used for the density conservation equation, we get: 



D 
Dt 



B+0 + + B-Q- 



e- 


.dS = 


f dS ± 












Dt 





Tn-u s B ± 5 s ± n.u ± B ± 8 



± H±XS 



(65) 



Now, taking into account that the terms proportional to 6 s , and V<5 S must vanish independantly from each 
other we will write the Faraday law for the bulk and for the interface as: 

• terms proportional to 9^: 



DB d 
~Dt 



= -Vx# 



(66) 



Let us stress that the total time derivative appears quite naturally because the continuous bulk media 
are moving 



• terms proportional to 5 s : 



-j^- r 2HB s u s .n +B s V s .u s + B+n.{u+ - u s ) - &-ft.{it - u s ) = -V x E s - n x (E+ - E~) (67) 

new term 

There appears now a drift flux of the inductions B^ , which did not show up in [Til ], since they considered 
stationary bulks. 

• terms proportional to \75 s : 

n x E s = (68) 
The equation for the Ampere law, written in the local form is: 

3D -* 

VxH= — +j (69) 
at 

However, we have again to go over to the integral form and consider thus instead of ([69]) 33Ll34|: 

(V x H)dS = — J D.dS + J j.dS (70) 
Applying the usual method to the Ampere law, we obtain now the following equations: 



terms proportional to 9^: 



terms proportional to 5 s : 



VxiJ^^+f (71) 



— n x (B+ - §-) + —V x B s -fix (M + - M~) - V x M s = e —^- 
D S P 



+ D s \7 s .u s + D + n.(u + - u s ) - D~n.{u- - u s ) + f - 2HD s u s .n (72) 



now term 



n x M s - —n x B s = (73) 
Mo 



terms proportional to W5 S : 



The two first Maxwell's equations (|55|) and (I59|) . describing the surface interaction with the electro-magnetic 
field are exactly the same as in [lg. But for the two others (|€>T[) and (|72^) . we note the presence of a new term, 
proportional to the mean curvature and of a drift of induction and displacement. At this stage, we do not 
know how big these terms are. All we can say is that they are theoretically there. 

Furthermore, Gogosov and coworkers of 0] consider the situation of a ferrofluid without electric field, thus 
they consider only the second |5^|) and the last Maxwell equations (ITTj) deduced respectively from ((57]) and 
in (|69|) where they neglect any electrical sources. They used a development in a power series of a small 
parameter, chosen to be the ratio of the width of the interface to a characteristic length measured along 
that interface which is supposed thus to be of a very large extent. Using such a perturbation method, they 
established the analog of the usual boundary condition for the induction and the magnetic field, but in the 
absence of a drift term, since they assumed continuity of the normal velocities along the surface. Gogosov's 
results [5;] show a dependence on the mean curvature to appear in the boundary conditions, as a consequence 
of the deformation and the dynamics of the interface. 



3. Continuity equation for the charge density conservation 



In its integral form, the continuity equation for the charge density is: 

^ J ( Pe )dV T = - J V.jdVr (74) 

Using the decomposition |T]) for the current j and the charge p e with Vr — S V + (B V~ as before, and using 
the properties © and |(HJ), we obtain the following results: 



±. 



• terms proportional to 9 

^ + ^ = (75) 

• terms proportional to <5 S : 

+ ptV s .u s - 2H{u.n) Pe + V.f + n.{j+ - J") + i o+n.(u+ - u s ) - p~n.(u- - u s ) = (76) 

• terms proportional to \7S S : 

n.f = (77) 

IV. EQUATIONS FOR THE ANGULAR MOMENTUM 
A. Total angular momentum 

Following [U H^, H(| , we define the total angular momentum of our system as : 

pK = rx (pu+p em ) + pd = pL + pCl + f x p em (78) 

We have already introduced the Poynting vector p em = D x B. The external angular momentum L due to 
fluid vorticity is defined as L = f x u, where r is the position vector and f2 is the internal angular momentum 
(spin density) introduced by individual particle rotation[lll[I3,[28j]. From [27j or [ll[ (derived by Kovac in the 
absence of electro- magnetic fields), we know that the dynamical equation for the total angular momentum 
takes the following integral form: 

J pKdV = J V.(f x f)dV + J V.YdU + J pCdV + J fx pF ex dV (79) 

where pC is the external torque applied by the system environment. The tensor Y is called stress couple 
tensor by Rosensweig [HI and convective flux of angular momentum by Kovac This tensor Y and the 
tensor T, using (j46|) . are defined as follows (see for more details [13,[2^| ): 

Y = r)'(VQ + (Vf>) T ) + 2£'Z.(V x Q) + A'(V.f2)I 

(80) 

f = f + T em + up em = T + 2£Z.(lu - O) + T em + up em 

Where rl is the shear coefficient of spin viscosity, A' is the volume coefficient of spin viscosity and £' is the 
vortex spin vorticity. We will apply the usual decomposition |T]) for the physical quantities pK, p em , Q, T, Y, 
respectively. We can rewrite (f79"|) as follow (using Vr — V + V" S); 

' dp±R±Q± ■ v.(K± P ±u±e±)) dv + f( W*' sa )) ds + [ Uk^ DdS 



dt ' v ' I J \ Dt I J V ) Dt 

J V.(rxT6 ± )dU + J V.(f x TS s )dV + J V.(Y ± e ± )dV ± + J p ± 6 ± £dV ± + J f x p ± ± F ex dV ± + 

J V.(Y s S s )dV + J P s S s CdS + J fxp s 5 s F ex dS (81) 



This equation becomes: 



/ y Ik 6 +^-(^ ± /° ± " ±e± )j + / \^^ + ^-(^P S K s )-2p s HK s {u s .n)j 5 s dS = 
J V.{rxf@ ± )dV + J V.(r x TS s )dV + J \7 .{Y ± <d ± )dV ± + J p ± Q ± £dV ± + J r x p ± Q ± F ex dV ± + 

J V.(Y s 5 s )dV + J P s 6 s CdS + J rxp s 5 s F ex dS (82) 

For the right side, we use the properties © and §S§ to obtain the following results: 

• bulk phases: 

dp±K± + V.Q^i?^) = V.(f x f ± + Y ± ) + p ± r x F ex + p ± C (83) 
at 

This is exactly the equation for the angular momentum used by Rosensweig [13, HH 

• interface: 
dp s K 



dt 



-V \(p s K s u s )-2H{u s .n)p s K s = n.u s (p+ K + -p~ K~)-n.(p + K+u+ -p~ K~u-)+V .(rxT +Y S ) 



+ n.(rx(T -T )) +n.{Y + -Y~) + p s £ + p s rx F ex (84) 

• boundary conditions (terms proportional toV<5 s ): 

n.(r x (f S + T s em + u s f em ) + Y s ) = (85) 

The two last equations extend the one found by Kovac [ll[ to a deformable interface (which might be flat) 
and in conjunction with an electro-magnetic field 

B. External angular momentum 

The equations for the external angular momentum L = r x u for each phase are obtained directly from the 
corresponding equations (|42|) and ([43|) as a starting point. Taking the vectorial product of (|42|) and ([43|) by r, 
the resulting equation is: 

^ + V.G^f^) = fx V.(T ± + Tfj -rx^+rx (p±F ex ) (86) 
But, using the property [i~7l ]: 

r x (V.T) = -V.(T x r) + Z : T (87) 

we obtain the following: 

^ + VVi^) = -V.(T ± x r + T± m x r) + Z : (T ± + T±J -fx^+fx (p±F ex ) (88) 
For the interface we start from the equation: 

() '' 1 '-V .{p s L s u s )-rx(2p s H(u s .n)u s ) = f x (V.T S )+p- L~ n.{u- -u s )- p+ L + n.{ii + -u s )+rx (n.(T + -T~)) + 



0f 



f x (n.(T+ m - T-J) +fx (V.T^J + n.u s (r x ft m -fx ftj -fx^+fx (p°F ex ) (89) 



Using (|5T[) and the relation 
one obtains finally: 



dp s L s 
dt 



'x (n.T) = -(n.T) x 



V.(p s L s u s ) -2H(u s .n)p s L s 



(90) 



p-L~n.(u- -u s ) - p + L+n.(u+ - u s ) - V.(T x r ) + Z : T — (n.(T — T ) x 



V.(T^ m x r) + Z : T* m - (n.(T+ m - T^)) 
together with the following normal boundary condition: 



3'D S —* 
xf- n.u s {f^ m - p; m ) xf-fx + f x (p s F e: 



dt 



x n.(T + XJ + u s .n(f x f em ) = 



(91) 



(92) 



C. Internal angular momentum 

Applying the B.A.M theoretical derivation, Kovac (see [ll|) considers a system including angular momentum 
effects, with external as well as internal angular momentum. But he does not take into account any electro- 
magnetic effects. Thus in his work (Tl| . the total angular momentum reduces to: K = fi + L, which simplifies 
the obtention of the equations for the internal angular momentum: he substracts the equation (|83p and (|84[) 
for the external angular momentum, from the corresponding equation (|88[) and (|9ip for K. 

When electro- magnetic fields are present, the definition of the total angular momentum is different as we 
recalled just above (see ((75|) (27l. I2H |30| . Then, we follow an approach analogous to the one of R.Rosensweig in 
[I?} combining it with the BAM derivation to obtain an equation describing the internal angular momentum 
balance in the bulk phase and in the surface. From the definition (|78p for the total angular momentum, (|79p 
becomes: 



fxpu + fxp em + pn dV = V.{fxT)dV + V.YdV+ pCdV + fxpF ex dV (93) 



D 
Dt 

We consider individually the different terms appearing in the left side of (|93p . The first term is: 



~Dt ^ ' '"'"'^ ' 



P ± §- t (rx(u±S ± )) 



dV 



-^(fxp s u s 5 s 



dS+ [fx p s u s 5 s 



DdS 
Dt 



(94) 



where we have by definition u — -j^ , 



— I f x pudV = 
Dt 1 



r x p 



w so that u x D 
Du ± S ± 



Dt 



dV + fx 



■jyf — 0. We obtain the following expression: 
dp 3 u s S s 



dt 



+ V.(p s u s S s ) - 2H(u s .n)p s u s S s 



dS (95) 



Now, applying the usual BAM derivation (see |T])up to ([2"9"|) and (J3D])), the second term in the LHS of ([9"3"|) 

gives: 



£ J f X p ern dv = ^- t J fx (^ m e±)dv + ^- t J fx (f em s s )ds 

Using systematically the decomposition (|57|. last expression can be rewritten as follows: 
D 



(96) 



Dt 



'x p em dV = 



z : (^e*) + f x ( -^g— + v.(u ± ^ m e ± ; 



dV 



dV-2 (Hu s J s {fxf em ))dS 



(97) 



The decomposition ((57)) applies also to the right side of (f§3"|) , beginning with the first term, rewritten as: 

V.(r x f )dV = [ \f X (V.T) - Z : fl dV (98) 



We have 



V.(f x T)dV = I r x V.(f ± 9 ± ) - Z : (f ± 6 ± ) + r x V.(f V) - Z : (TV 



If we put all the results in equation (|9"3"|) . we can rewrite it as follows: 

fx (p ± ^(^e ± )-v.(TV) + ^- P ± te ± l '/r 



(99) 



/ X [ dpS ^ 5S + v.(p s u s 6 s ) - 2H(u s .n)p s u s 5 s - V.(T>) - 2Hu s n (f em )S s + - p s F e ^ 



V 9* 



eft 



Z:(T 9*) + V.(Y ± 6 ± ) + V.(Y S J S ) +p ± £6 



p s £-Z:T 5 s 



(100) 



The RHS of last equation is equal to zero since it is the vectorial product of r by the momentum balance (|43p 
and (|42|).for the bulk and for the interface. We get the following: 



^ f [ p ±n±Q±)dv + ^ lunr^uis 



Z:(T Q ± ) + V.(Y ± Q ± )+V.{Y s S s )+p ± Ce d 



dV+ 



p s C — Z : T 6 s 



dS 



(101) 



This is the equation for the internal momentum, which we were looking for. Then, the usual properties and 
procedure leads us to describe each bulk by: 



V.[/fn ir) = V.Y - Z : (T -fT^J+p^C 



(102) 



The equation at the interface has the form: 

dp s & 
dt ' 



V.{p s n s u s ) -2H(u s .n)p s n s 



= V.Y S + n.(Y + - Y") - Z : T - Z : T* m - /9 + 0+n.(u + - u s ) + p-fl-R.(iT - u s ) + p s C (103) 
and the boundary conditions are 

n.Y s = (104) 
V. CONSTITUTIVE EQUATIONS FOR THE RELAXATION PROCESSES 



To close the system of equations, we need to have dynamical equations for the magnetization, M and for the 
electric polarization P. For example, in the case of a ferrofluid, one can consider it as a homogeneous one 
component stable solution [13, However for such a system oscillating far from its stationary state, the 
collinearity assumption between the field H and the magnetization M does not hold and the equation of state 
of the form M = M(T, H) = \ H (where x 1S a constant), is not adequate anymore [13 [H, [2t| [3(| • It becomes 
mandatory to have a dynamical equation for both quantities M and P in order to describe the relaxation 
process which occurs in this kind of situation. Although a lot of works has been developed to obtain these 
equations, mostly in relationship to the molecular origin of the relaxation process [13, HI], there is no real 
agreement on what is exactly the right equation to use at the macroscopic level. This is why, we will consider 
the most common phenomenological equations for both polarization and magnetization relaxation processes, 
appearing in the literature. 



A. Relaxation equation for the magnetization 



For the modeling of the relaxation process of the magnetization, we consider the Shliomis relaxation relation- 
ship, [28j |. generalized for compressible fluids [2§|: 

^ + V.(mM) = fi X M + —(M ~ M) (105) 
at T m 

where r m is the magnetic relaxation time constant and Mq is the value of the magnetization at the equilibrium 
given by the usual Langevin equation [T3|. We assume that the value of these two quantities can be different 
for each phase of the system. Using the usual BAM derivation, we obtain then: 



• For the bulk: 

^— - + V.(u ± M ± ) = Q ± x M ± + —^(Mq - M±) (106) 



We have considered the generalized form of the equation just for mathematical easiness. Nevertheless, 
in this work, we consider the bulk phase as a incompressible medium, then the equation (|106|) can be 
simplified as follow: 

nw± i _ 

— — = fi± x M ± + — (Mq - M±) (107) 

Dt T m 



For the interface: 



— — + V.(u a M s ) - 2H(u s .n)M s + M + n.(u+ - u s ) - M~n.(vr - u s ) = s x M s + — (Af s - M s ) (108) 

We notice that we do not have an equation proportional to V<5 S . It means that the boundary conditions for the 
magnetization are completely determined by the boundary conditions obtained from the Maxwell equations. 



B. Relaxation equation for the polarization 

For the polarization we can use a similar equation as the one for the magnetic relaxation process (see also 

My- 

^- + W.(HP) = -(Po-P) (109) 
dt r e 

Where r e is the electric dipolar relaxation time and Pq is the value of polarization would be in equilibrium. 
The usual procedure leads to the following set of equations: 



For the bulk: 



°— + V.(^P±) = ^(P ± - P±) (110) 



As explained for the magnetization equation, for the bulk phase we can have the following simplified 
form: 



Dr - P 1 ) I 111 I 



Dt 



• For the interface: 



BP S - ^1^^ 

— + V.(u s P s ) + P + ri.(u + - u s ) - p-n.(u- - u s ) - 2H{u s .n)P s = — (P S - P s ) (112) 

As for the magnetization, there is no term leading to supplementary boundary conditions for the polarization. 
However, let us stress that in usual magnetic fluids the electric polarization is small, and in practice the 
constitutive equation for the polarization can be reduced to: 

P » Pq = Z e e E (113) 

where £ e is the electric susceptibility. To obtain the corresponding equations for the bulk phase and for the 
interface is, in this case, straightforward. 



VI. EQUATIONS FOR THE ENERGIES OF THE SYSTEM 



In the last section wc derived the dynamical equations for the physical variables of our system: the momentum, 
pu, the electromagnetic field, B, D, the angular momentum, pK. Now we want to obtain the corresponding 
equations for the different kind of energies which exist in our system: we will begin with the kinetic energy 
and the electromagnetic energy, after that we will consider the equation for the angular momentum energy, 
and we will finish with the conservation equation for the total and the energy the internal of the system. We 
have to notice that if we want to do this work we need to reformulate the dynamical equations, using the two 
following relations (p?Tj) and (J3U), that we reproduce here: 



dp A 
dt 



+ V.(pFu ± ) = 



(114) 



do 8 

+ V.(p s M s ) + p + n.(u+ - u s ) - p-n.(u- - u 8 ) - 2H(n.u 8 )p s = 

dt 



(115) 



With these relations, we will be able to rewrite the equations in a lagrangian form, as we show below. Let's 
consider the following expression: 



dp s F s dF 
P • V.{p s F s u s ) - 2H(rt.u s )p s F = p s ' 



dt 



dt 



rr.vr ) f- ( ^ + v.{ P s u s ) - 2H{n.u s ) P s ) j i i :i c) 



This relation hold for any quantity F, be it a scalar F, a vector F or a tensor F. Now, if we take into account 
the relation (|115[) , we obtain from (|116[) : 



Q p s F s 

dt 



V.(p s F s u s ) - 2H(n.u s )p s F 



,DF S 
Dt 



F s (p + n.(u + - u s ) - p n.(u 



')) 



(117) 



We have a direct relation between the eulerian derivative, and the lagrangian derivative we note that 
in the last case, the curvature term doesn't appear anymore. We have a similar result for the equation inside 
the bulk: 



dp s F ± - . . ^. . ± DF ± 
1 \7.(p ± F ± u ± ) = p ± 



Ot 



Dt 



(118) 



In the next sections we will use (|117j) and (|118|) to transform the equation in a lagrangian form, which will 
simplify the calculation for the equations for the different kind of energy. 



A. Equations for the kinetic energy 

With (fTTT)) and l[TTS]l . equations (gSJ) and can be rewritten as follow: 
• Equations for the bulks: 



^-(DxB) ± +p ± F ex 



(119) 



Equations for the interface: 



Dt 



= p-(n.(vr -u s ))(u- -u s ) - p + (n.{u + - u s )){u+ - u s ) + V.f +n.(f -f ) + p s F e , 



(D x B)+ - (D x BY 



(120) 



• boundary conditions: 
These conditions remain the same as before: 



n.(f +T s em ) + u s .n(DxBY = 



(121) 



The equation for the kinetic energy, defined as U C i n = ^u 2 , can be obtained by taking the scalar product of 
the last two equations with the vector u. 
We have for the bulks: 

P^-u* = (V.(t ± + Tt m ).u±) - + P ± F ex .u^ (122) 

We use use the relation: 

(V.T).u = -V.(Tu) + T : W (123) 

to obtain the formulation: 

P± Wt^ ] = + O^) + (** + Tf m ) : W± - + p±F ex .u± (124) 

At the interface, the equations have the form: 

(V.(T S + T s e j) .u s + (n.(T + - T")) .u s + (n.(T+ m - T em )) .u s + (n.u s )& m - P^ m ).u s - ^ff -U s (125) 
which become, using the same relation (|167p : 

ps Wt (y) = p ~ K* (ir _ ^' ~ -" s " p+ " " S)) ^ ~ ^ MS ~ 

V. ((f S + T s em )u s ) + (f S + T em ) : Vtf + (n.(T + - T - )) + (ri.(T+ , - T em )) 

(n.u s )(# TO - jCJ.tf* ~ ^ff (126) 

For the next part, it is useful to rewrite the equations in a more conventional form, in term of dissipation 
functions and fluxes. 



• For the bulk, we have: 



at 

Where we defined the current: 



+ V.^Z^J + V. J± + <4 = (127) 



^ n = (T ± +T^).tf t (128) 



and the dissipation function: 



4, = ^ff - P ± ^-« ± - (* ± + O : W± (129) 



• For the interface, we obtain 
dp s U; 



at 

with the current 

and the dissipation function 



V.(p s u s W cm ) + V.J s cm + o s cm - 2H(u s .n)p s W cin = (130) 
fs in = (T s +Tl m )M s (131) 



<in = ^.u s -p s F ex .u°~(f S + Tl m ) : W s -(n.(f + + T+ m -(t +T e " m ))).^ -(n.u s M m -fe m ).u s + 

U^ n n.{p-(u- - u s ) - p+(u+ - u s )) + u s . [p + (ti.(u+ - u s ))u+ - p-(n.{vT - u s ))u~] (132) 

The procedure for the next part is the same as here, then we will show directly the results, skipping the details 
of the calculation. 



B. Equations for the electromagnetic energy 



To obtain the equations for the conservation of the magnetic energy for the bulks and the interface, we will 
strat from the equations (|6*6"1 ITTj) for the bulks and from equations (|67|72|) for the interface. 

• Equation for the bulks: 
We multiply equation ([6"6"|) by and equation (fTTj) by we obtain 



DB 



± 



= -^.(V x E±) 



Dt 



(133) 



^.(V x H ± ) = E± 



DD i 
Dt 



(134) 



Substracting equation (|134p from equation (|133|) and using the identity V.(D x H) = H.(V x D) 
E.(V x i/) we obtain after some manipulations 



D(H ± .B ± ) D0 ± .E ± ) 



Dt 



V.^ x H z 



DH± DD* 
-i -D = —^—.B +D ± 



Dt ' " v ' ' J Dt ' 

Using the phenomenological relations D = /iq(H + M) and D = CpD + P, we can obtain 



D 

Dt 



■ V-C-E* x i? ± ) + j^.D* 



DP d 
Dt 



■.D = 



DAP 



Dt 



= 



This equation can be rewritten in an analog form as for the kinetic energy: 



dU. 



■± 

cm 



dt 



V.{^Ut m )+V.Jt m + ai m = Q 



with the following definition for the electromagnetic energy in the bulk 
and with the dissipation function 



. - DP ± p., DM± 



Dt 



Dt 



and with the current 



Jt m =E ± xH ± 



(135) 



(136) 



(137) 



(138) 



(139) 



(140) 



• Equation for the interface: 
We start from the equations (fBTf and (|72| . We follow exactly the same method we used for the bulk 



phases and we obtain the following equation: 



Dt 



D P s -> - \ I D M s 
E s . | -4r— - 2H(u s .n)P s + P s V s .u s + H s .[ -|- 2U{u s .ft)M s + M S V S 



Dt 



Dt 



H s .{nx(E+ -E-))-E s .(nx{H+ ~H-))+n.{u+ -u s )(H s .B+ +E S .D+)-n.(u~ ~u s ){H s .B~ +E S .D~) = 

(141) 

= (142) 



Which can be rewritten in the following way 



dt 



V.{u s W em ) + V.J| m + a\ m - 2H{u s .n)U { 



with 



(143) 



(144) 



a% m = -2H{u s .n)W em + (V s .u s )W em + E s .f + E s . - 2H(u s .n)P s + P s V s .u s J + 

(D M s -, -, \ ^ -> - -. - - 

— ^- 2H{u s .n)M s + M s \7 s .u s + H s .(n x (E+ - BT)) - E s .(ft x (H+ - H~)) 

+ n.{u + - u s )(H s .B+ + E S .D+) - n.(tT - u s )(H s .B~ + E S .D-) (145) 

• boundary condition: we have the simple following condition 

n.J s em = (146) 

In the definitions (|139[) and (|145[) of the dissipation functions appear the total time derivative of the polarization 
and of the magnetization. In order to obtain a closed set of equations for the electromagnetic energy itself, 
we can use the constitutive equations for the polarization and for the magnetization. 

C. Equations for the rotational kinetic energy 



To obtain this equation we follow the same procedure as we did to obtain the equation for the kinetic energy. 
We define the internal energy for the rotation as (U ro t = h& 2 )- The calculation leads to the following results: 
For the bulk 

P± \^»T = (^- Y± )-^ - (Z : (** + Tfj) .ft* + p ± C.n ± (147) 



For the interface, we find the following expression 



p 



2 Dt 



{V.Y s ).n s - (z : (T 3 + T s em )\ .& + P s C.tt s + (n.(Y + - Y - )) 



((&- - n s ).&^j (p-n.{u- - u s )) - ((&+ - fi s ).d s ) (p+n.{u+ - u s )) (148) 

In term of dissipation function and flux, these equations become 

• For the bulk: 

+ V.^tt^) + V. Jtt + °tot = (149) 

with the current: 

J± i =Y ± .0 ± (150) 

with the dissipation function 

att = (Z = (T ± + Tf m )) .fr - Y ± : Vft± - p ± C.VL ± (151) 

• For the interface, we have: 

dp s U, 



1 < V.(p s W rot u s ) + V..P rot + a s rot - 2H{n.u s )p s U s rot = (152) 



at 

with the current 



J? ot =Y 8 .fr (153) 
and the dissipation function has the following form 

o* rot = (Z : (T S + Tjj) - Y s : VO s - p s £.& - (n.(Y+ - Y")) 
fr. \fl-p-n.(u- - u s ) - n + p+n.(u+ - u s )] + U s rot n.{u- - u s ) - U s rot p + n.{u+ - u s ) (154) 



VII. CONSERVATION OF THE TOTAL ENERGY AND INTERNAL ENERGY 



We start with the total energy, p£, which is defined as 

p£ = pU in + pUdn + pU rot + U em (155) 
This quantity is conserved, therefore the conservation equation has the form: 

^- + y.(p£u) + y.J e = (156) 
at 

with the current flux, J e : 

J e Jq -\- J cin ~l~ Jem ~l~ J rot (157) 

Now, we can rewrite it as: 

j p£dV T = J V.J e dV T (158) 

Applying a similar method as we did for the equation for the mass conservation we have the following results: 
for the bulks we have the well known equations: 

^—+V.{ P ± u ± £ ± )+V.Jt = Q (159) 

This equation is just the equation (|156|) . with the bulk's variables. For the interface, the equation takes the 
following form: 

+ V.J e s + £+n.{u+ - u s ) - £-n.(u- - u s ) + n.(J+ - J~) - 2H(n.u s )p s £ s + V.(p s u s £ s ) = (160) 
with the boundary condition: 

J e .n = (161) 

To obtain the equation for the internal energy, pUi n , we need to start from the equation for the total energy 
and substract from it the expressions of the equations for the other quantities, which leads to the following 
expression: 

+ V.(piEY cin ) + V. J q = a em + a cin + a rot (162) 

The use of the usual decomposition procedure leads to the following results: 

• For the bulks: 

^™ + V.G^O + V. J± = af m + a±„ + af ot (163) 

• For the interface: 
?^+V.(p s i! s U s c J+VJs+nM s (p-Ur^^ 

(164) 

• Boundary condition for the current 

ft. J s q = (165) 

VIII. ENTROPY PRODUCTION AND FLUXES 

In this last section we consider the entropy production of our system. The first step will be the derivation of 
the equation for the evolution of the entropy and in a second step we will compute the different fluxes. As we 
did before we can write the equation for the entropy, S in a integral form: 



D 
151 



pSdV + / V. J s dV = / <£ s dV (166) 



Applying the usual method, we obtain the following set of equations: 



• For the bulk: 



dp ± S ± 
dt 



+ V.(jf+p ± u ± S ± ) = $ 



(167) 



• For the interface 
dp s S 



dt 



V.(J S S + p s u s S s ) - 2U{u s .n)p s S s = $ s s + p-S-ft.{u- - u s ) - p + S + n.{u + - u s ) (168) 



• boundary condition: we have 



ft.Jt = 



(169) 



We have to specify the entropy currents Jf- and J| and also the dissipation functions & s s and We start 
with the well known relation between the internal energy and the entropy: 
For the bulk: 



DS ± p ± DUf n p ± Dp ± E ± DP ± H± DM ± 



+ 



Dt T± Dt p ± T ± Dt T±' Dt 



Mo 



T± ' Dt 



(170) 



For the interface: 



P 



,D S S S P s D s W ln P s Hp' , E s I) J'" , //" D S M S 



Dt T s Dt p s T s Dt T s ' Dt ' Mo ~T S ' Dt 



Using the definition of the lagrangian derivative, we have: 
For the bulk: 



(171) 



^+V.(p±u±Ut) 



p± Dp± t £± DP ± , i?± DM ± 

Mo — —RT- ( 172 ) 



p±T± Dt T± ' Dt ^ T± ■ Dt 



For the interface: 



dt 



\7.(p s u s S s ) - S s p-n.{u- - u s ) + S s p+n.{u+ - u s ) = 



1 



dt 



+ V.(p s u s W m ) - 2H(n.u s )p s U? n - U s lnP -n.(U- - u s ) +W inP + n.{u^ - u s ) 



p- D s p s E s D S P S IP D S M S 
p s T s Dt + T 7 ' Dt Dt 



Using the equations for the internal energy we obtain 
For the bulk: 



(173) 



dp ± S ± 
dt 



V.ip^S*) = 



1 



-V. Jf + af m + af m + a± ot 



p ± Dp ± S ± DP ± , itrA . 



p±T± Dt T±' Dt 



Mo 



H ± DM± 
T±' Dt 



For the interface: 



dp s S s 
dt 



+ V.(p s u s S s ) - S s p-n.{u- - u s ) + S s p+n.{u+ - u s ) = 



1 



V.J* - n.u s (p-Ur n p+U+) n.(p + u + U+ p~iTUr) 



+ 



[< m + < ln +<ot-'2n{n.u s )p s U s m -U s m p n.(u - u s ) +W m p + n.{u+ - u s )} 



p s D s p s E s D S P S 



H s D,M S 



p s T s Dt T s Dt 
Now, we can take into account the following relation: 



Mo 



T s Dt 



(175) 



Lv.J q = V.(^) + 7 Lj q .VT 



then we can obtain the equation: 
For the bulk: 



dp ± S ± 
dt 



J* 



V.(p ± u ± 5 ± ) + V.(^) = T 



-^±(VT).J±+a± Tl 



(176) 



cin ' rot 



p ± Dp ± 
p ± T ± Dt 



DP ± H ± DM 



± 



T± ' Dt ^ u T± ' Dt 
comparison with (|167p leads to the following identification 



Ate-^- (177) 



Jt = 4: (178) 



The left side of (|177j) can be seen as the dissipation function $ ± 



For the interface: 



+ V.(p s u s S s ) - S s p-n.(u- - u s ) + S s p+n.(u+ - u s ) + V.(^) = 



- — (VT). J« - H.u s (p-Ur - p+U+) - n.(p + u+U+ - p-u-Ur n ) 



+ ^ [-2H{u s n) P s Uf n + u s em + u s cin + a s rot - Uf nP n.{u - u s ) + W mP + n.{u+ - u s )] 

p s D s p s E s D S P S H s D S M S , . 

S J—-\ ._f 1- n ._f 179) 

we can compare this expression with (|168l) , and we obtain 

Jt = §i (180) 

The left side is the function $^ + 2Hp s S s (u s .n). In the expressions we have obtained here, some terms which 
are function of the mean curvature appear, there are absent from the previous work for the same reason as 
explained before. If we compare our results with the results in we notice some differences coming from 
the fact that in [5], the drift terms, as we said before, have been neglected. Thus, the expressions we have 
obtained here are a generalization of the previous results. 




IX. CONCLUSION 



In this paper, we have shown how the effect of the time dependent curvature can be taken into account in the 
framework of the theory built by D. Bcdeaux, A.M. Albano and their coworkers. We obtained very general 
equations for a system of two immiscible liquid phases with an interface, for which system, the presence of 
singular electromagnetic field is taken into account. In the limit of the simplifications which were done in 
previous works, we recover the results of [l[ — Q. We have also stressed some features of the equations in their 
lagrangian form for which the non linear term disappears leading to simpler but complete equations. We end 
up with a unified theory for the mathematical treatment of such systems, since various previous approaches 
have been shown to lead to identical results. In future works, we plan to develop the Onsager relations and to 
show what is exactly the effect of this non linear term. Indeed, most authors have not given any illustration 
of the effect of this term or explained how it should be or not important. 
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